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Abstract 

We use realistic two- and three-nucleon interactions in a hybrid chiral-perturbation- 
theory calculation of the charge-symmetry-breaking reaction dd — ► an to show that 
a cross section of the experimentally measured size can be obtained using LO and 
NNLO pion-production operators. This result supports the validity of our power 
counting scheme and demonstrates the necessity of using an accurate treatment of 
ISI and FSI. 



1. For most purposes, hadronic isospin states can be considered as charge symmetric, i.e., 
invariant under a rotation by 180° around the 2-axis in isospin space. Charge symmetry (CS) is 
a subset of the general isospin symmetry, charge independence (CI), which requires invariance 
under any rotation in isospin space. In quantum chromodynamics (QCD), CS implies that 
dynamics are invariant under the exchange of the up and down quarks [1]. However, since 
the up and down quarks do have different masses (m u ^ ma) [2,3], the QCD Lagrangian is 
not charge symmetric. This symmetry violation is called charge symmetry breaking (CSB). 
The different electromagnetic interactions of the up and down quarks break CI. Observing the 
effects of CSB interactions therefore provides a probe of m u and m^. 

Two exciting recent observations of CSB in experiments involving the production of neutral 
pions stimulate our attention. Many years of effort led to the observation of CSB in np — > dn 



Preprint submitted to Elsevier Science 



9 February 2008 



at TRIUMF. The CSB forward-backward asymmetry of the differential cross section was found 
to be An, = [17.2 ± 8(stat) ± 5.5(sys)] x 10~ 4 [4]. In addition, the final experiment at the IUCF 
Cooler ring reported a very convincing dd — > an signal near threshold (a = 12.7 ± 2.2 pb at 
T d = 228.5 MeV and 15.1 ±3.1 pb at 231.8 MeV) [5]. The dd -> an reaction violates CS 
since the deuterons and the a-particle are self-conjugate under the CS operator, with a positive 
eigenvalue, while the neutral pion wave function changes sign. 

The study of CSB 7r° production reactions presents an exciting new opportunity to learn about 
the influence of quark masses in nuclear physics, and to use effective field theory (EFT) to 
improve our understanding of how QCD works [6]. This is because chiral symmetry of QCD 
determines the form of pionic interactions. Electromagnetic CSB is typically of the same order 
of magnitude as the strong one, and also can be handled using EFT. 

The EFT for the Standard Model at momenta comparable to the pion mass, Q ~ m 7T , is chiral 
perturbation theory (xPT) [7]. This EFT has been extended [8,9,13,11,10,12] to momenta 
relevant to pion production, Q ~ y/m n M with M the nucleon mass. (For a review and further 
references, see Ref. [13].) 

EFT with the operators of Ref. [14] was used to correctly predict the sign of the forward- 
backward asymmetry in np — > dn° [9]. For the dd — > air reaction, we surveyed various mecha- 
nisms using initial-state plane-wave functions and simplified final-state wave functions [10]. In 
this simplified model, we found that the formally leading-order (LO) production mechanism 
is suppressed through symmetries in the wave functions and studied other mechanisms. The 
contributions from next-to-next-to-leading-order (NNLO) diagrams are too small to account 
for the observed cross section — a cross section of only 0.9 pb was found. We also included 
short-range pion emission, which contributes at N 4 LO through contact vertices whose strengths 
are a priori unknown. We used resonance saturation, by means of CSB effects in Z-diagrams, 
as motivated by a successful phenomenological model [15] of the charge-symmetry-conserving 
(CSC) reaction pp — * ppn°. For the simplified wave functions, we then found a cross section of 
the observed order of magnitude. 

Our aim here is to take advantage of recent significant advances in four-body theory [16,17] 
that allow us to include the effects of deuteron-deuteron interactions in the initial state, and 
to use bound-state wave functions with realistic two- and three-nucleon interactions. The cal- 
culations presented here are hybrid: the pion-pro duct ion operators are constructed using EFT, 
but the nuclear interactions used to obtain the wave functions are not. No calculation of this 
kind can be considered to be completely well-founded unless the operators and wave functions 
are constructed from the same convergent EFT. However, the interactions do include one-pion 
exchange, so their long-range behavior is founded in EFT. Moreover, we employ several poten- 
tials to gauge the sensitivity of the various production operators to the shorter-ranged parts of 
the interaction. 

The present study does not include all diagrams appearing at NNLO. A complete analysis 
demands the inclusion of loop diagrams. Their evaluation requires a careful treatment of di- 
vergences as pointed out in Ref. [11], and understood only recently [12]. For technical reasons, 
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(a) (b) (c) 

Fig. 1. Diagrams of np — > d-rr ; the solid circle indicates CSB. 

photon exchange is so far only considered in the final-state wave function. In this letter, we 
concentrate on the important effects of initial- and final-state interactions (ISI and FSI), and 
anticipate that the use of the present incomplete set of operators should be sufficient to get 
order-of-magnitude estimates and to demonstrate the technique. 



2. We summarize the power counting of the previous study [10], which contains explicit ex- 
pressions for the operators. At LO, there is only one contribution, represented by Fig. la: 
pion rescattering in which the CSB occurs through the seagull pion-nucleon terms linked to 
the nucleon-mass splitting. This contribution stems from the chiral transformation properties 
of the quark operators that generate CSB, which are two: i) the up-down mass difference, 
which breaks chiral symmetry as a component of a chiral four- vector; ii) electromagnetic quark 
interactions, which break chiral symmetry as components of a chiral anti-symmetric rank- 
two tensor. In lowest order, there exist two seagull operators involving a nucleon interacting 
with two pions, one of which is neutral. Their strengths are determined by the quark-mass 
and electromagnetic contributions to the nucleon mass splitting, 5m N and 5m N , respectively. 
SrriN is proportional to e{m u + m d ), where e = (m u — md)/(m u + m d ) ps 1/3, while 5m N is 
the fine-structure constant a times a typical hadronic mass. The only existing constraint on 
these two terms is 5m^ + Sm^ = M n — M p = 1.29 MeV and the model dependent estimate 
Stun = — (0.76±0.3) MeV based on the Cottingham sum rule [21]. Verifying the theory requires 
that the two terms be constrained independently. The most natural reaction to study is irN 
scattering. Ref. [2] predicted a significant difference between the 7T°p and ir n scattering lengths 
that is not presently observable, as discussed in Ref. [18]. Effects of these terms in the nuclear 
potential [19] are relatively small or suffer from other unknowns as in 7id scattering [20]. This 
leaves the investigation of CSB in the two reactions, np — > dn and dd — > an , as very promising 
possibilities. For definiteness, in this paper we use the central value of the estimate from the 
Cottingham sum rule. The leading diagram, Fig. la, is 0[eml/(f%MQ)], where f n = 92 A MeV 
denotes the pion decay constant and Q ~ y/m^M a typical momentum. 

We refer to this contribution as "pion exchange". For completeness, we distinguish the parts 
proportional to Sm^ and Sijin and denote the contributions by M.pe = -MpE,Sm N + -M.pE,Sm N - 

There is no NLO contribution. At NNLO, suppressed by 0(m^/M), there exists a recoil cor- 
rection of the LO term (labeled Ai rec = M. r ec,&m N + M-recSm N )- Its strength is also determined 
by 5rriN and 5m^. Therefore, its contribution allows us to estimate the size of NNLO contribu- 
tions. The recoil correction to the ttNN vertex is linear in the energy of the virtual pion. Such 
operators were studied in Ref. [22] for the reactions AW — > NNir, and we use the prescription 
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provided there and applied in Ref. [10]. Demonstrating the validity of this recipe for a four-body 
environment deserves further study. 

At the same order new parameters appear. In particular, a term arises in which a one-body CSB 
operator (oc Pi + /3 3 ) is sandwiched between initial- and final-state wave functions, as illustrated 
in, e.g., Fig. lb. We refer to this as the one-body term (A^i&). The terms (5\ = 0(em^/M 2 ) 
and /6?3 = 0(a/ir) arise from, respectively, the quark-mass-difference and electromagnetic con- 
tributions to the isospin-violating pion-nucleon coupling. Neither f3\ nor (3% can be extracted 
from experiment yet. To allow us to provide numerical results we estimate these terms by 
modeling [23] @ 1 by n-r] mixing, see Fig. lb, 

^ = g v {AH\v)/K (1) 



where (ir°\H\ri) = —4200 MeV 2 is the vr-r^-mixing matrix element [24], and g v = g V NNfir/M 
the 77- nucleon coupling constant. An early analysis [25] using one-boson-exchange potentials 
in NN scattering gave g 2 NN /4iT = 3.86 (used in Ref. [9]), but the data show little sensitivity 
to r] exchange and high-accuracy fits can be achieved [26] using g 2 NN /4ii = 0. Indeed, the 
possibility of a vanishing coupling constant had been raised earlier. The detailed analysis of 
NN total cross sections and pp data using dispersion relations [27] found that g 2 NN /An = 0. 
This is consistent with extractions from the nucleon pole in the amplitude ttN — > r]N that give 
[28] 0.5 > g 2 NN /4n > 0. Photoproduction reactions on a nucleon [29] (see their Fig. 2) yield the 
small value g 2 NN /4ir = 0.1. To be consistent with our earlier study [10], we use g 2 NN /4ir = 0.51 
for the results shown below, but also examine the effects of using g 2 NN /4n = 0.10. Both values 
are roughly consistent with the size expected using power counting arguments [23]. We assume 
the sign predicted by SU(3) symmetry, as in Ref. [10]. 

The effects of electromagnetic interactions as well as strong CSB were included in computing 
the ct-particle wave functions, where the former effect is dominant. These interactions generate a 
small isospin T = 1 component of the wave function that enables a non-zero contribution of CSC 
production operators. To estimate the effects of the admixtures, we calculate the production 
matrix element using the CSC counterpart of diagram Fig. lb (referred to as M. WF ). 

A number of other CSB mechanisms enter at N 3 LO or higher, including additional loop dia- 
grams and short-range interactions. The lowest order where four-nucleon contact interactions 
start to contribute is N 4 LO, that is, 0(m 1T /M) below NNLO. To estimate their strength, 
Ref. [10] evaluated certain tree-level contributions as indicated by Fig. lc, which represents 
the exchange of heavy mesons (a, ou, p) via a Z-graph mechanism, with ir-r] mixing generating 
CSB at pion emission (A4 a , M. u and M. p ). Another Z-graph (labeled as M. pu} ) arises in which 
the CSB occurs in the heavy-meson exchange via p-uj mixing along with strong pion emission 
at the vertex. The Z-graphs are believed to be important because their inclusion leads to a 
quantitative description of the total cross section for the reaction pp — > ppn near threshold 
[15]. Our present results use the coupling constants and parameters of Ref. [10], see their Ta- 
ble I. However, in the future it will be necessary to reassess the procedure in light of recent 
developments concerning the treatment of divergences in EFT loop diagrams [12]. 
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3. The various mechanisms generate pion-production kernels that are sandwiched between 
final- and initial-state wave functions to provide a transition matrix element M.. We restrict 
our analysis to T d = 228.5 MeV, as the effects of a small change in energy are captured mainly 
by the change in the phase-space factor. The cross section is related to the matrix elements by 



a = 4.303 pb 



M 10 4 fm 2 



(2) 



We present our new results in stages. First we introduce realistic bound-state wave functions, 
while continuing to use the plane- wave approximation (PWA). The techniques to solve the four- 
body problem have been presented by Nogga et al. [16]. To be specific, we present results using 
both the AV18 [30] and CD-Bonn 2000 [26] two-nucleon potentials combined with a properly 
adjusted Tucson-Melbourne (TM99) [31] three-nucleon force. The combination guarantees that 
the a-particle binding energy is reproduced with high accuracy. Additional calculations using 
the Urbana-IX [32] three-nucleon potential resulted in essentially identical results and will be 
presented elsewhere [33]. 

Table 1 summarizes our results for the transition amplitudes that add to M., labeled according 
to the various mechanisms described above. The one-body term is predicted rather model 
independently. Using these matrix elements for the one-body operator leads to a cross section 
of 10 — 13 pb, which is accidentally in good agreement with the experiment. Compared to our 
toy- model calculation [10], we find an increase of the cross section by a factor of 10, showing 
that the high- momentum tail of the wave function is important, especially for the one-body 
term. Using the smaller, but also realistic, coupling g^ NN /An = 0.10 would reduce the resulting 
cross section by a factor of 5. The one-body term is formally subleading. However, the toy- 
model calculation showed that the pion-exchange term is suppressed due to the symmetry of 
the a-particle wave function. This result persists for the realistic a-particle wave functions: the 
amplitude does not vanish exactly, but remains smaller than the one-body term. This term is 
quite sensitive to the chosen nuclear interaction, pointing to sensitivity to the short-range part 
of the potential and to the small components of the a-particle wave function. 

Since the LO term is suppressed and the one-body term is not well constrained by resonance 
saturation, it is interesting to look at the pion-recoil term. Its parameters are better determined 
than Pi, since they are related to the nucleon mass difference and the Cottingham sum rule. 
Our calculation may therefore give a trustful estimate of the size of the NNLO contributions. 
The results are a rather model-independent amplitude of approximately 1/3 the size of the 
one-body term and are in line with the power counting. 

In contrast, we find that all the Z-graphs give unexpectedly large contributions, especially the 
p-uj exchange operator. Also, the contributions add constructively, so that their sum tends to 
overwhelm the one-body term. This model of resonance saturation thus gives results in vast 
disagreement with the power counting. 

CSB effects on the final-state wave functions (AAwf) are smaller than pion recoil and insen- 
sitive to the chosen nuclear interaction, indicating that these terms are well constrained using 
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Table 1 

Complex dd — > air amplitudes at = 228.5 MeV in units of 10~ 4 fm~ 2 . PWA denotes the plane-wave 
approximation. ISI results also include the initial-state interaction. 





PWA 

CDB+TM99 AV18+TM99 


ISI 

CDB+TM99 AV18+TM99 


A/I r> t? 

A4 dip i 

Mp E 


0.35 -0.07 
0.06 -0.01 
0.41 -0.08 


-1.51 + i 1.87 -0.76 + i 0.74 
-0.28 + i 0.35 -0.14 + i 0.14 
-1.79 + i 2.22 -0.90 + i 0.88 


J L rec 


0.41 0.34 
0.08 0.06 
0.49 0.40 


-0.81 + i 0.74 -0.63 + i 0.59 
-0.15 + i 0.14 -0.12 + i 0.11 
-0.96 + i 0.88 -0.75 + i 0.70 


M lb 


1.76 1.60 


-2.51 + i 1.84 -1.94 + i 1.60 


M a 
M u 
M p 


0.46 0.31 
0.51 0.38 
0.24 0.15 
1 17 0.87 


-0.56 + i 0.64 -0.32 + i 0.42 
-0.53 + i 0.44 -0.35 + i 0.34 
-0.33 + i 0.34 -0.18 + i 0.19 
-1.32 + i 1.51 -0.84 + i 1.07 


M W F 


-0.15 -0.14 


+0.51 - i 0.13 +0.41 - i 0.14 



phenomenological interactions. 



4. The next step is to present the effects of including the ISI. The correct treatment of this 
involves the solution of the four-body scattering problem at center-of-mass energies greater 
than 200 MeV. In spite of the tremendous progress achieved in recent years on the solution of 
the four-nucleon problem [34,17], advances in obtaining exact solutions are limited to energies 
below the four-particle breakup threshold. To understand our pion-production reaction it is 
necessary to go beyond the distortions obtained through an effective optical-model potential 
fitted to the elastic dd scattering data. This is because important pion production occurs in 
which the deuterons interact, break up, and then emit a pion. At very high energies the use 
of Glauber approximation is justified, but in the threshold energy regime for pion production 
the wave length associated with the relative d + d on-shell momentum is close to the size of 
the deuteron. Therefore we obtain an approximate solution of the Yakubovsky [35] equation 
for the four-nucleon scattering wave function that is made up of two terms: the first involves 
the bound-state wave functions of the two deuterons times a plane wave describing the relative 
motion between them; the second requires the breakup of one of the deuterons followed by the 
three-body scattering of the N + d system into the three-particle continuum in the presence of 
the remainder spectator nucleon. 

Such an approximation is based on the lowest-order terms in the Neumann series expansion of 
the four-particle Yakubovsky equation, leading to the following expression for the scattering 
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wave function, 



|*Po) ~ l^iJJC, * G 0g 5 PP J«#>> , (3) 

3 ip 

where p denotes one of the seven two-body partitions, four of (3) + 1 type and three of (2) + (2) 
type, and i is a pair interaction that is internal to p; j is both internal to p and p . The 
initial-state wave function component \<pj°) carries the appropriate bound-state wave function 
components of the target and projectile times a relative plane wave between their respective 
center of mass. As usual, p specifies the two-body entrance channel, 5 PPq = 1 — 5 PPq , G is the 
four- free-particle Green's function and t { the t-matrix for pair i. If p corresponds to a 2 + 2 
initial state, then p can only be a (3) + 1 two-body partition and Ufj is the solution of the three- 
body Alt, Grassberger and Sandhas (AGS) [36] equation for the three-particles that make up 
subsystem p. The first term in Eq. (3) corresponds to the initial-state wave function; the second 
term requires the breakup of one of the bound pairs followed by the scattering of either one of 
the particles from the remaining bound pair, leading to four free particles in the continuum. 
The particle-pair scattering into the continuum takes place in the presence of the fourth one, 
and therefore, by energy conservation, the total energy available is the total four-body center- 
of-mass energy minus the relative kinetic energy of the fourth particle relative to the center of 
mass of the other three. Thus the four-body scattering wave function we construct contains all 
orders in the pair interaction but also three-particle correlations in first-order perturbation at 
all possible energies that are consistent with four-particle energy conservation. 

For four identical nucleons Eq. (3) may be written as 



l*£> = \M + ^= (1 + P - PuP + P)(l - P 3 4)|^ 2 ' 3)4 ), (4) 

where P = P12P23 + P13P23 and P = P13P24 are permutation operators whose appropriate 
combination generates the 6 (12) components of the Yakubovsky wave function that are of 
2 + 2 (1 + 3) type. The symmetrized d + d initial-state wave function [10] \4>dd) is given by 



\(f>dd) = ^= (1 + P ~ P34P + P) 6(12) 6(34) Exp(12 - 34), (5) 

where £d{ij) is the deuteron wave function for the pair (ij), and Exp(12 — 34) represents the 
relative plane wave between the two deuteron pairs. The second term in Eq. (4) mandates the 
use of a specific choice of wave-function component for nucleons 1, 2, 3 and 4, 



l^ 2 ' 3)4 ) = G (12,3|[/oWI(12)3)a34), (6) 
where connectivity increases from left to right and z = E — ^ k 2 + iO, k being the relative 



momentum between nucleon 4 and the center of mass of (123). The breakup operator U = tG U 
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and the corresponding matrix element (12, 3 1 Uq(z)\(12)3) represents the scattering of nucleon 
3 from the bound state of (12) leading to three free nucleons in the continuum where nucleons 
1 and 2 are last to interact through their respective t-matrix t. The operator U is the AGS 
three-body scattering operator that satisfies the integral equation 

U = PGu 1 + PtG U, (7) 

from which one calculates Nd elastic scattering amplitudes. The permutation operator P is the 
same as used in Eq. (4) and corresponds to the sum of the two cyclic permutations of particles 
1, 2 and 3. We extract from Eq. (4) the 3 P partial wave in the entrance channel to compute 
the threshold cross section for dd — > an . 

Details regarding the numerical solution of the scattering problem and evaluation of the pion- 
production matrix-element will be presented elsewhere [33]. Here, in Table 1, we simply present 
the results for the transition amplitudes, which acquire imaginary parts due to the presence of 
the initial-state interaction. Generally, we observe a significant enhancement of all contributions 
to the amplitude. 

It is immediately apparent that the pion-exchange term, which is supposed to be LO, is now of 
the size of the NNLO terms considered here, namely the pion-recoil and one-body operators. 
It still shows a sizable model dependence, which could visibly influence our final result for the 
cross section. A more consistent treatment of nuclear interactions and production operators 
will be necessary in future. 

The pion-recoil and one-body terms remain relatively model independent. Both, therefore, 
can serve as an order-of-magnitude estimate of the cross section. Our results for these matrix 
elements correspond to cross sections between 4.5 and 42 pb. Again, the one-body contribution 
is larger than the pion-recoil term. Both would come close to each other for the smaller choices of 
the strength of the one-body term. Our explicit calculation shows that the NNLO contribution 
provides a strength consistent with the experiment. We stress that the strong enhancement due 
to initial-state interactions and higher-momentum components of the a-particle wave functions 
are necessary to find NNLO contributions of the required size. 

The ISI also enhances short-range contributions from Z-graphs, but by far-smaller amounts. 
We still find relative contributions much larger than expected from the power counting and, 
again, all the contributions add up constructively. One possible explanation would be that there 
is simply no convergent EFT for the reaction dd — > an and the Z-graphs would need to be 
included as done in this work. However, if this was true the value of the empirical cross section 
would be the result of subtle cancellations amongst various terms from very different origins, a 
very unlikely coincidence. The more likely interpretation is that the power counting works, but 
the Z-graphs simply provide the wrong model to estimate the four-nucleon operators. If this 
is the case we may even drop them all together from our investigations, since they are of high 
order. 
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5. This paper extends our earlier study of the reaction dd — > air by using realistic wave func- 
tions for the four-nucleon ISI and FSI. We also provide numerical estimates for some diagrams 
that lead to a cross section of the right order of magnitude supporting the power counting given 
the suppression of the LO. In addition, the present results allow us to identify a few issues that 
deserve further study (in addition to the inclusion of all diagrams up to NNLO). 

• Given the dramatic influence of initial-state interactions, it is of paramount importance that 
new experimental constraints be obtained for the deuteron-deuteron interactions in the energy 
region close to the pion-production threshold. Besides data on elastic dd scattering, also data 
on other pion-production reactions with the same initial 3 Po state are needed. The most 
obvious examples are the CSC reactions dd ^ 3 H/HeiV7r recently measured at COSY [37], 
which are an interesting alternative to elastic dd scattering because only a few partial waves 
contribute in the entrance channel. 

• Another important issue is to better understand the role of the Z-graphs used to estimate 
the size of four-nucleon operators. As shown above, their size is much larger than indicated 
by their N 4 LO power-counting order. Thus, it is necessary to reassess the procedure in light 
of recent developments in EFT, especially concerning the treatment of divergences in loop 
diagrams [12]. Given the large number of experimental data, especially for pp — * ppn°, much 
insight can be obtained. 

• Finally, we found some dependence on the chosen nuclear interaction model. This is prob- 
ably related to the presumed inconsistency between the chosen nuclear interaction and the 
production operators. Such an inconsistency can only be resolved by applying nuclear inter- 
actions based on chiral perturbation theory [7,38] that need to be extended to higher cutoffs 
as outlined in Ref. [39]. 

The central result of this letter is that the inclusion of ISI and FSI enhances the contribution of 
NNLO diagrams of the production operator enough to able to account for the measured cross 
section. Together with the insight that the LO is suppressed, this supports the EFT approach to 
this reaction. But it also demonstrates that a careful treatment of the nuclear effects is required 
for a final analysis. 

In view of a planned measurement of the reaction dd — » air at higher energies at COSY [40] 
— where p waves will be relevant — and of the experimental determination of Aj b [4], we will 
have an increased database on cross sections, which will allow us to disentangle the various 
contributions. Although much remains to be done before any precise statements about the 
values of the parameters Sm^, Sm^ can be made, we are now at threshold of understanding 
how the light-quark masses makes a difference in nuclear physics. 
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